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Abstract. For the Fermi-Pasta-Ulam chain, an effective Hamiltonian is con- 
structed, describing the motion of approximate, weakly localized discrete breath- 
ers traveling along the chain. The velocity of these moving and localized vibra- 
tions can be estimated analytically as the group velocity of the corresponding 
wave packet. The Peierls-Nabarro barrier is estimated for strongly localized 
discrete breathers. 

1. Introduction. Discrete breatliers (DB) hiave tlie defining properties of being 
spatially localized and having time-periodic dynamics. They are also called in- 
trinsically localized, in distinction to Anderson localization triggered by disorder. A 
necessary condition for their existence is the nonlinearity of the equations of motion 
of the system, and the existence of discrete breathers has been proved rigorously 
for some classes of systems ^1 01 01 ^1 ^2 ^] . In contrast to their analogs 
in continuous systems, the existence of discrete breathers is a generic phenomenon, 
which accounts for considerable interest in these objects from a physical point of 
view in the last decade. In fact, recent experiments could demonstrate the existence 
of discrete breathers in various real systems such as low-dimensional crystals 
antiferromagnetic materials |221, Josephson junction arrays jH], molecular crystals 
|7] , coupled optical waveguides ) and micromcchanical cantilever arrays [23 • For 
a review on the topic see 

More than a decade ago, discrete breathers have been observed in numerical 
simulations of, among others, the Fermi-Pasta-Ulam (FPU) system (5311201: a chain 
of nonlinearly coupled masses. Although existence of DB as exact periodic solutions 
was proved already some years ago for a large class of oscillator networks ^1 
(not containing the FPU system), a proof of their existence in the FPU system has 
been obtained only recently ^1 Ej . 

A generalization of the concept of DB are so-called traveling DB, spatially lo- 
calized solutions which travel along the chain. At least as approximate solutions 
they have been observed numerically in the FPU system. A proof of their existence, 
however, has not been achieved so far, and in fact it is doubted that traveling DB 
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exist as exact solutions in this system. The theoretical analysis of approximate 
traveling DB is still an open problem which will be tackled in this article. We con- 
sider traveling DB as long-living transient structures which can be described by an 
effective Hamiltonian as presented in [I7| . We distinguish three types of localized 
structures, characterized by their degree of localization, which can be found in the 
FPU dynamics. Examples of these different types are plotted in Figure 1, where 
from type I to type III the degree of spatial localization is increasing. Type I is 
similar to a standing wave. It is not a proper DB, as its localization depends on the 
system size, but its energy envelope is sort of localized in a finite system. Type II 
is exponentially localized, but, as the width of the structure is not that small with 
respect to the intersite distance, it is called a weakly localized DB. Type III depicts 
a strongly localized DB, similar to those first discovered numerically [211. three 
of these structures are observable in FPU systems as (transient) traveling localized 
objects, called traveling discrete breathers, and we will distinguish in the following 
between these three types when deriving their effective dynamics. 
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Figure 1. Three types of traveling DB in FPU. Type I: weakly lo- 
calized standing wave. Type II: weakly localized DB. Type III: strongly 
localized DB. The red line in the upper plot represents the local energy. 



The outline of this article is as follows: Section |2] provides some basic definitions 
regarding the FPU model. In Section O the method of the effective Hamiltonian 
is briefly recalled, providing the necessary ingredients for its application. Section^ 
succinctly recalls some results from ^7jj describing the motion of DB of type I. In 
Section |S1 an effective Hamiltonian is constructed for a family of DB of type II, 
providing a good framework to understand the traveling motion of weakly localized 
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DB in FPU chains. Finally, in Sectional the case of strongly localized DB is treated. 
Here, instead of calculating the effective Hamiltonian explicitly, the so-called Peierls- 
Nabarro barrier is computed, which provides already useful information about the 
moving properties of DB of type III. 

2. FPU model. Let us consider a chain of oscillators whose dynamics is deter- 
mined by the Hamiltonian 



(1) 



with momenta p„ g ]R and positions a;„ € R. V and W characterize, respectively, 
the on-site potential and the (nearest-neighbor) interaction potential. The FPU 
model is defined by the choices 

V^O and W{x) = ^x'' + ^x'^ + ^x^, (2) 

2 3 4 

where fc, a, /? G R. Note, however, that the theory of the next section applies to a 
more general class of models. 

3. Method of the effective Hamiltonian. The method of the effective Hamil- 
tonian has been exposed in | 17| in a general framework. It is briefly recalled here 
before applying it to the FPU model in the subsequent section. 
Let us define a traveling DB of as a function of the form 

Xn{t) = u{A,n — ct,Cjt), (3) 

where the function u{A, r, s) is increasing in A, and such that u(0, •, •) ~ 0. It is 
localized in r (e.g., |u(-,r, •)! < e"'^'''' for some A > 0) and periodic in s, u(-, -,5-1- 
27r) = - ,5). To be slightly more general, one may think of ct being replaced by 
a function Q{t) monotonic in i, where Q = c + f{t) and / is a function of small 
amplitude, periodic with a frequency which is small compared to uj. The point is 
that, in a moving frame with velocity Q, one sees a discrete breather at rest with 
frequency uj. 

It is generally believed that solutions of the form Q do not exist as exact solu- 
tions in generic systems, although a proof of this conjecture has been accomplished 
only for a certain class of systems 5 . However, it may be possible to get approxi- 
mate traveling DB of the form Q when the following conditions can be met: 

1. Find a two-parameter family of approximate (non-traveling) discrete breathers 
of (^3) of the form 

Xn{t) = u{A,n- Q,u;t), (4) 

where A characterizes the amplitude of the DB and Q its "position" . Typ- 
ically, for a chain having some translation symmetry, and for fixed A, there 
are two such positions where there is an exact DB: a "site-centered" DB with 
integer Q, and a "bond-centered" DB with half- integer Q. By interpolating 
between these exact solutions, approximate DB can be constructed for any 
real Q. 

Next, consider the family of functions 

xi'^^^^'"') =u{A,n-Q,ujt- kn). (5) 

The idea is that there may be a "Doppler effect" of the traveling DB, oj — uj = 
kc, which can be described by adding a phase shift kn to tot, where k is called 
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the momentum of the DB. The point is that ^ stiU describes approximate 
DB for small enough k. 

Compute, for each member of the family x[^'^'^\ the symplectic area 



a 



piA,Q,U).iA.Q.U) (6) 



with period T = lixju. In the case where p = x, like for the Hamiltonian 
one can write 

a = 2T{E^,^{A,k,Q)) , 

where (•) denotes the mean value over one period and iJ^in is the kinetic 
energy. This calculation should result in a symplectic area a proportional to 
A, or related to A by a simple functional relation. If this is not the case, 
the parametrization in jSjl should be changed, indexing the family in terms 
of {a,Q,k). The fact that a is a proper choice as a parameter is explained 
in | 15II17[05| . In short, this can be justified by stating that a is an adiabatic 
invariant of the dynamics • 

3. Compute the symplectic form dp„ A da;„ restricted to the coordinates 
(fc, Q). This is achieved by evaluating 

^'^Q = t], \^^-dQ~-dQ-dk)'^'' 

which should be non-zero for any (k,Q), and the restricted symplectic form 
reads TkQ dk A dQ. 

4. Then, an effective Hamiltonian is computed as the mean (over one period) of 
the original Hamiltonian along the family of DB, 

Hf (fc,Q) 1^ H[{xt'^ '^\ptQ '^^})dt, (8) 
and, finally, the dynamics of {k, Q) is given by 

TkQ dQ 

Q ^ ^'-^^ (10) 

TkQ ok 

For small enough fc, and assuming dkK^ik, Q) = for fc = (otherwise shift 
k k — ko to have this) one may expand the effective Hamiltonian so that it 
takes on the form 

providing that m(Q) be non-zero and finite. The latter defines an effective 
inertia, although m(Q) may be negative. In fact, as discussed in ^Ij (especially 
in examples of Section 6 therein), the sign of m{Q) can be interpreted as the 
sign of an effective anharmonicity of the breather solutions. The second term 
in is a periodic function (of period 1) which defines a Peierls-Naharro 
potential for discrete breathers, as introduced and discussed in El El- 
In these papers it is shown that the critical points of V°^(Q) correspond to 
exact discrete breathers of the full model. In the simplest situation there are 
two types of critical points associated, respectively, with the stable discrete 
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breathers and the unstable ones. The difference of at these points, say 
Q'"-* and Q^^\ is called the Peierls-Nabarro barrier 

AEp^ = |V'=*(Q(")) - V'=*f(Q(''))|. (12) 

The absolute value in this equation is necessary in case of negative m{Q). 
For m[Q) positive, the minima are associated with stable DB, and the max- 
ima with unstable DB, otherwise the stability is just reversed. In any case, 
the smaller this barrier is, the higher is the mobility of a traveling discrete 
breather, and the better it is in general described by the effective Hamiltonian. 

In the next sections we apply this scheme to the analysis of the mobility of various 
types of localized solutions in the FPU model. 

4. Traveling localized v^raves (type I). In this section, a particular type of 
localized solutions is considered, termed type I in the introduction, which have the 
shape of a standing wave. They do not constitute proper DB, as the localization 
depends on the system size. The method of the effective Hamiltonian was applied to 
this type of solutions in |17j . Here, we summarize these results in order to compare 
them to those for the proper DB solutions obtained in the following sections. 

Consider an FPU chain with, say, an odd number M = 2iV-|- 1, iV € N, of masses 
and periodic boundary conditions. There exists, in the linear regime (a = /3 = 0), 
a family of standing waves 

xJt) = -^=\fA (-1)" sin f — 71 - Q] cos(wt) (13) 

whose envelope is localized in space, as it is a sine function with wavelength equal 
to twice the chain length (see Figure^ type I). Moreover, it is observed numerically 
that the localization strength is enforced when the nonlinearity is turned on, a, /3 ^ 
0. This family (|13|l complies with as it is parametrized by A and Q, while 
uj = 2 sin (37-) is fixed. Such a wave can be put into slow motion when a and /3 
are non-zero. 

To "add momentum" to this standing wave, it is convenient not to follow 0, 
but rather to consider the expression 
1 



Xn{t) = —j= ^/ A + k sin [tot — cm + Q) — ^ A — k sin {ojt + an — Q) , (14) 

where a — 2N+1 ^^'^ ^ taken as the momentum parameter conjugate to Q. 
The motivation for favoring this parametrization is related to the structure of the 
linearized FPU chain and has been explained in |17| . 

With these settings, the symplectic area ^ can be computed, yielding a = 
A/2tt. The computation of t^q from Eq. Q gives t^q — 1, and so the effective 
symplectic form corresponding to the coordinates (fc, Q) is simply canonical. Finally, 
one computes the averaged Hamiltonian JHJ, obtaining 

Hf ik) ^ CO A +f^P^{3A'-P), (15) 

from which the equations of motion 

3 uj^ 

are deduced. Therefore, to this order of approximation, there is no Peierls-Nabarro 
barrier, and the motion of a traveling localized wave \Vi\ is a slow drift at constant 
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speed MQ/ir} This result is confirmed by direct numerical simulation of the FPU 
system in the weakly nonlinear regime, using H14|) as initial condition and plotting 
the variable Q versus time t (Figure The prediction of Hl()|) turns out to be valid 
for small times t. The discrepancy encountered for larger t might be reduced by 
considering higher order terms in H15|l . 



(a) (b) 




10 20 30 40 50 200 400 600 800 1000 

time time 



Figure 2. Time evolution of the coordinate Q. Straight lines 
indicate the predictions of the effective dynamics (|16|l . The pa- 
rameters of the FPU model in ^ are chosen as a = 0, /3 = 1, 
A = 0.05, Qo = — f ; and feg = —0.01. Panel (a) is an enlargement 
of a part of panel (b), corresponding to the time lapse [0, 50]. 



5. Weakly localized discrete breathers (type II). The first proof of the ex- 
istence of DB in FPU type systems is due to Livi, Spicci, and MacKay ^1] for a 
diatomic chain of alternating masses m and M with large mass ratio m/M . Exis- 
tence of DB in monoatomic FPU systems [as defined in Q and |2Jl] has been proved 
only recently, either using a variational approach or elaborating a technique of 
central manifold reduction |10[ [TT] . The latter existence result obtained by James 
will be used as a starting point in order to construct a family of DB of the form 
Q. In his formulation |nil^2i James works with the force variable 

Un = W{X„ - Xn-l) (17) 

which is one-to-one related to the displacements a;„, provided that W is invertible. 
Suppose 

W"iO) = ^ = 1. (18) 

Then the main result of TIT can be formulated as follows: In the FPU system, for 
u! — utQ > sufficiently small, there exist DB of the form 

y„ ~ (-l)"u„ cos(tjt)-hC'(t^-Wo), (19) 
where O denotes the Landau order symbol, and m„ satisfies the recurrence relation 

Un+l + Un-1 - 2u„ = K^Un - "fU^ + h.O.t. (20) 

with the condition |m„| — > as n — )■ ±00. Higher oder terms (h.o.t.) have been 
neglected in (|20ll . Equivalently, u„ (with rt e Z) must be a homoclinic orbit of 



Note that the speed depends on the system size. 
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MO = „„...r.:.. „M (23) 



the two-dimensional map (Mn,M„-i) (m„+i,m„) defined in (|20(l . In this map, the 
parameters are defined as 

k'^u^-A and ^ = (21) 

where 7 has to be positive in order to guarantee existence of DB. So the assumption 
Lo"^ ^ ojq — i imphes that <C 1 in In addition, let us suppose that ^ 

7|ii„p. Then, can be approximated by the differential equation 

v" ^v- (22) 

with Un = -^v{nK). Equation (|22|l has a family of solutions homoclinic to u = 0, 
namely 

V2 

cosh[K(^ — q)] 

parametrized by g G R. Hence, one deduces an approximate analytic expression 
for ifTOI) . reading 

/2" cos(ijjt) 

ynit) ~ i-irnJ ; ^ 24 

V 7 cosh[K(n — (y)J 

which represents a family of approximate DB of the form As an example, 

a member of this family is plotted in Figure ^ (type II) for k = 0.2. Finally, 
substituting ujt by cut — kn in the cosine function with fc w tt, provides the family 
of functions 

[2 coslujt — kn) . . 

y„ t ~ kJ 25 

V 7 cosh[K(n — Q)\ 

of the form which can serve as a good starting point for the construction of an 
effective Hamiltonian for (approximate) traveling DB. Equation H24|) is recovered 
for k = TT. The family H25|l is indexed by three independent parameters, Q G R, 
|fc - 7r| < 1, and k < 1. 

5.1. Change of coordinates. The Hamiltonian ^ is written in canonical coordi- 
nates (a;„,p„). It turns out to be convenient to work in relative coordinates (p„, J„) 
defined by 

Pn = Jn-Jn+1, (26) 

and in the following wc will rewrite the FPU Hamiltonian with y = 0] as well 
as the family of DB H25I) in these coordinates. We assume the original coordinates 
to satisfy the conditions 



< CXD, 



^\Xn-Xn-l\ < OO, 

n 

hm x„ = 0. (27) 

n — ^ — OO 

For the new variables, these relations imply 

E|J„-J„+i| < OO, 

n 

J2\pn\ < OO. (28) 
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In addition, it will be assumed that lim„^_oo Jn — 0. Then the change of vari- 
ables l|26f) is invertible. For example, for given {(pn, </n)}neZj one can compute 

n 

Pn^Jyi-Jn+l and Xn= Rk- (29) 

A;— — oo 

Moreover, the transformation is canonical since 

dpn A dXn = ^ (dJn A dXn-1 + dJn A dpn - dJn+1 A dXn) 
n n 

= ^dJ„Adp„. (30) 

n 

In the new coordinates, the FPU Hamiltonian takes on the form 

H = J2liJn-Jn+lf + W{pn), (31) 

n 

and the canonical equations read 

Pn = "^Jn — Jn+1 ~ Jn-1, (32) 

Jn = -W'ipn), (33) 

or, equivalently, 

Pn - W^'(Pn+l) + W'{pn-l) - 2W'(pn). (34) 

In order to compute an effective Hamiltonian from H31I) we need to express the 
family of DB (|25|l in terms of the variables (p„, J„). This is achieved by considering 
the simple approximation 

Pn = {W')-\yn)^yn + 0{yl), (35) 

which is valid due to j/„ being small. The corresponding momentum J„ can be 
deduced from yielding 

Jn{t) = - I M^'(p„(s))ds-f J„(0) 
t 



(36) 



- / yn{s)ds -t- J„(0) 
Jo 

K sin(cijt — kn) 



w V 7 cosh[K(n — Q)] 



where J„(0) = —-\ - — f'^"^^"jQ^1 has been chosen. Note that assuming (p„, J„) 



cj y 7 cosh[K(n— Q)] 

to be small does not imply that Ansatz (|25() can be considered as a wave-packet 
obeying the linearized equations of motion. In fact, the linear limit 7 — > renders 
singular. 

5.2. EfTective Hamiltonian for weakly localized DB. In this section, we show 
that the effective dynamics of approximate traveling DB (|25|l can be described by 
the Hamiltonian 

H°*^(fc,g) = i-(3-cos/fc)-fC'(K2). (37) 
27 ' 
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This result is obtained by following the procedure recalled in Section |3 First, 
compute the symplectic area 



a = y JnPn^ dt. 



(38) 

Using H32f) . this expression can be written in terms of J„ only, reading 

a - 2 / V ( J,^ - JnJn+i) dt. (39) 

•^0 n 

Substituting H3t)|) in this equation yields 

AtTk'^ -yy ( 1 cosfc 1 



^^1 „ I cosh^[K(n - Q)] cosh[K(n + 1- Q)] cosh[K(n - Q)] 
This series can be estimated making use of the Poisson summation formula 

5]/(n-Q)=5]F(m)e2-™« (41) 

n rn 

where the coefficients F{m) of the Fourier series are the integer samphng of the 
Fourier transform 



F{v) = / J{x)e'^''"'^dx (42) 

J —OQ 

The advantage of this formula is that it is sufficient to keep only a few terms in the 
Fourier series ()41|1 when F{m) decreases rapidly for m — > ±oo. Here one obtains 

V ^ = - + '^'^ , 2. cos(27rg) + O [cxp (- — \\ (43) 

^cosh'Kn-Q)] K K^sm\i{^) ^ ^' I ^\ -J\ ^ ^ 

and 

y I = . (44) 

^ cosh[K(?i + 1 — Q)] cosh[K(n — Q)] sinhK 

For small k, which corresponds to our assumption of w — c^o being small, the second 
term in H43|l is exponentially small compared to the first one. Neglecting this small 
term and using H21|) , the area reads 



7(4+^2)3/2 

This result indicates that Ansatz H25|l does not define a family of DB whose area 
a is independent of k. Consequently, following Section |31 one should in princi- 
ple reparametrize this family by considering the change of variables (k, k, Q) — )■ 
(a, fc, Q). On the other hand, and for the sake of simplicity, one may assume k ^ n 
and, again, k small, to obtain 

a^2TT'^+0{K^,\TT-k\^) (46) 

at lowest order. In this case, a is proportional to k and no change of parametrization 
is needed. 

Next, following Section |21 the restriction tj-q dfc A dQ of the symplectic form to 
the coordinates (fc, Q), is determined. From {Tj), using a similar technique as for the 
computation of a, one finds at lowest order 

nQ = - + 0{K'). (47) 

7 
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Then, the effective Hamiltonian 

H'"" = ^ E - JnJn+1 + W {pn)] dt 

= ^(3-cosfc) + 0(K2)+V'=«(g) (48) 
27 

is computed. The effective potentiai defines a Peierls-Nabarro potential for 
traveUng DB I17| . It is a periodic function of period 1 whose first Fourier term 
reads Vb cos(27r(5), with 

7sinh(^)' 

and the so-called Peierls-Nabarro barrier is defined by 2Vq. As discussed above for 
Eq. H43|l . for values of k consistent with our hypothesis <C 1, such a term is 
exponentially small and therefore negligible. 

The velocity Q of the DB can be estimated from the effective equations of mo- 
tion 0, finding 

. 1 aH^ff sinfc , 

Let us compare this result with the group velocity of a wave packet with wavenumber 
k 

= cos - , (51) 



dfc V 2 



following from the dispersion relation of the linear chain uj(k) — 2 sin (-1) . Therefore, 
assuming |7r — fc| ^ 1, ifSUIl and are equal to first order in |7r — A;|, and one can 
write approximately 

In conclusion, the velocity of the approximate traveling DB given by (|25() is very 
close to the group velocity of the corresponding linear wave packet. Let us note that 
this observation was already reported in |22|, although in this reference Eq. (|52|l 
was not deduced, but rather taken as an assumption. The result is also confirmed 
by numerical integration of the equations of motion, illustrated in Figure |3Ja). 

Let us remark that Ansatz (|25|l does not behave like a linear wave packet. Firstly, 
it has already been mentioned that this Ansatz has a singular behavior in the linear 
limit 7—^0. Secondly, the maximum amplitude of a true linear wave packet would 
evolve in time like \pn\ oc For the weakly localized DB H25|l. however, the 

amplitude is approximately proportional to the area a, and is therefore predicted to 
be approximately constant in time. This difference in the time evolution is confirmed 
by numerical simulation. In Figure Efb) , the time evolution of the (inverse of the) 
DB amplitude l/max„ |p„P is plotted for two cases, namely with and without the 
presence of a nonlinearity, showing, respectively, a decreasing and a constant (in 
time) amplitude. 



6. Strongly localized discrete breathers (type III). In this section, we will 
not explicitly compute an effective Hamiltonian, although this should be possible. 
Instead, we will restrict ourselves to estimating the Peierls-Nabarro barrier directly 
from the properties of static (=non-travcling) DB. 
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Figure 3. (a) Numerical integration of the equations of motion 
of the FPU chain with initial condition H25() and k = 0.26, [3 — 0.5 
and k = TV — 0.4. A DB (solid) travels with a speed comparable 
to, but slightly larger than that of the corresponding linear wave 
packet (dashed), (b) The (inverse of the) amplitude of a trav- 
eling DB (solid) is oscillating around an approximately constant 
value, whereas the inverse of the amplitude of a linear wave packet 
(dashed) clearly increases linearly in time. 



6.1. Analytic estimate. Two kinds of strongly localized DB (called "type III" in 
Figure^ were found in the early numerical simulations of intrinsic localized modes 
of anharmonic lattices. Their shapes were approximately represented as 

< = ,0,-^,^,0,---) and = ,0,-i,l,-i,0,---), 

(53) 

known respectively as the Page mode (P), which is stable, and the Sievers-Takeno 
mode (ST), which is unstable. The approximations are used for large A^luq/P 
'13', and thus are valid in a regime different from the one for which Ansatz (|25|l 
holds. As an family of DB, interpolating between the P and the ST mode, we use a 
compact approximation (proposed by Kivshar 13 in a different context), namely 

Xn{t) = i-ir Au{n - Q)f{Lut), (54) 

where w is a truncated cosine defined by 

u(n) = I c°^(f") if l"l < 1' (55) 
[ otherwise, 

and / is a Jacobi elliptic function, e.g., f{ujt) — cn{ujt; k) (modulus n = \f2 in 13 ). 
Although it is known that the FPU system does not support compact solutions 
[H|, a convergence towards a compact function u [albeit not exactly towards the 
truncated cosine H55|l ] is observed in the large amplitude limit (= the limit of strong 
localization). Therefore, (|54|l is used here as a sketchy approximation from which 
the, likewise sketchy, ST and the P modes (|53|l are recovered for integer and half- 
integer values of Q, respectively. Equation H54|l defines a family of approximate DB, 
complying with the requirements of Q), and, following the procedure of Section |31 
it could be used to construct an effective Hamiltonian. Instead of doing so, we 
will proceed straight away to the estimation of the Peierls-Nabarro barrier from 
properties of the static DB, extending the results presented in |23| . 
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Let US first estimate the area of each member of the family H54|l . A simple 
substitution of this Ansatz in © yields 

aQ = 27Tu;QAl{imJ2^'in-Q)^ (56) 

n 

where Q — or Q — ^, and (•) denotes the mean value over one period. Inserting 
(I55|l . J2n u^{n — Q) = § for both the P and the ST mode. Therefore, the two modes 
have same area if 

LOp{Ap)Al = ujst{Ast)AIt, (57) 
(where the symbols ST and P represent, respectively, any integer or half- integer 
value Q). This equation is equivalent to saying that Ap is a function of ^st, e.g. 
Ap = A{Ast)^ and the function A can be evaluated at least numerically. Let us 
denote the corresponding periods Tst — 27r/a;sT(^ST) and Tp = 2tt /ujp{Ap). Then 
the Peierls-Nabarro barrier can be computed as the difference between the energies 

and 

^ ir'""^^''^^ = ^+A|(^c.g(/^)+Ap^a(/3)+A||/5(r)). (59) 

6.2. Numerical results. Numerical simulations were performed in order to check 
the validity of the approximations H58|) and (|59|) . and hence of the underlying as- 
sumptions (|54|1 and H55|l . As parameter values, we chose cjq = 1, a = 0, and (i — 1 
in the interaction potential |(2Jl. Discrete breathers can be calculated numerically up 
to machine precision, for example by means of a continuation of periodic orbits from 
the anticontinuous limit jl8| . We have computed four different quantities related 
to the mobility of a discrete breather, which are compared in Figure 0] 

1. Numerically exact Peierls-Nabarro barrier: Computing ST and P modes with 
the same value of the symplectic area a, the numerically exact Peierls-Nabarro 
barrier AiSpN = |i?p — ii^sxl can be obtained as the difference in energy between 
the two modes (dashed in Figure QJl. 

2. An approximation to the Peierls-Nabarro barrier: Equations (|58|l and (|59|l 
allow to approximate the Peierls-Nabarro barrier by computing the difference 
between their right hand sides (solid in Figure 0)). The values {p) = i, 
{P) = 0, and (/■*) = I have been used, which are found when / is assumed to 
be a cosine function. The values of the oscillation periods Tst and Tp, as well 
as the (maximum) DB amplitudes Ast and Ap , have been extracted from the 
numerical data. 

3. A second approximation to the Peierls-Nabarro barrier: Computing discrete 
breathers numerically, one notices that the oscillation periods Tst and Tp 
corresponding to the same value of the symplectic area a are clearly different, 
and Tst > Tp. In fact, it is observed that the main contribution to the 
above approximation (2.) of the Peierls-Nabarro barrier AE'pn stems from 
the kinetic part. 

This quantity is used as a second approximation of the Peierls-Nabarro barrier 
(dotted in Figure 0}. 
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4. An upper bound on the depinning energy: The aim of the computation of 
the Peierls-Nabarro barrier is to obtain an estimate of the energy necessary to 
"depin" a discrete breather (i.e., to kick a static DB such that it starts trav- 
ehng along the lattice). An upper bound on the depinning energy (points in 
Figure can be obtained from a simple numerical experiment which consists 
in kicking a static DB with a certain momentum distribution, make it evolve 
in time by means of numerical integration of the equations of motion, and 
see if it starts traveling along the lattice. This procedure, of course, depends 
strongly on the momentum distribution applied, and a rough numerical opti- 
mization scheme was employed in order to approximate the minimum energy 
necessary for depinning. 



10000 





Figure 4. Different approximations for the Peierls-Nabarro barrier 
ABpN versus the symplectic area a. Dashed: numerically exact Peierls- 
Nabarro barrier, solid: approximation (2.), dotted: "kinetic" approxi- 
mation (3.), points: upper bound on the depinning energy (see text for 
details) . 



Since Ansatz (|54ll is supposed to be valid in the limit of large A^loq / (3 with our 

choice of parameters), large values of the DB amplitude (and hence the symplectic 
area a) have been considered. However, the quality of the approximations does 
not change significantly with increasing a, and the relative error between exact 
result and approximation seems to be more or less constant. The quality of the 
approximations is reasonable, albeit not excellent, differing from the exact Peierls- 
Nabarro barrier roughly by a factor of two. 

When computing numerically the Peierls-Nabarro barrier A£'pN(a) — \Ep{a) — 
EsT{a)\ for some value of the symplectic area a, we were surprised to find Ep(a) > 
-Est (a). Typically, a stable mode is expected to correspond to the minima of the 
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Peierls-Nabarro potential, and an unstable mode to the maxima. Finding the reverse 
situation, we conclude that the motion of a traveling DB in the FPU system is 
effectively described by the dynamics of a particle with negative effective mass in 
the Peierls-Nabarro potential. 

The upper bound on the depinning energy we found in the numerical experiment 
differs from the Peierls-Nabarro barrier roughly by a factor of two. A more refined 
optimization of the momentum distribution employed in the depinning experiment 
might further improve the match. 

7. Conclusions. Effective Hamiltonians are constructed for various types of DB 
which can be found in the FPU model. These types, labelled type I, II, and III 
in Figure ^ are characterized by their degree of spatial localization. Type I is 
not a true discrete breather, as its localization depends on the system size, but 
it represents a localized wave which can be set into motion. The Peierls-Nabarro 
barrier is found to be zero in first approximation, and the localized wave travels at 
constant speed. An explicit estimate shows that this speed is not a group velocity 
and depends on the system size. DB of type II are called weakly localized. They 
are small in amplitude and can be analytically approximated by employing a recent 
result by James proving existence of DB in the FPU system. In this case, 

the Peierls-Nabarro barrier is found to exist, but it is negligible, being exponentially 
small with respect to the DB amplitude. DB of type II are found to travel at a 
speed approximately equal to the group velocity of the corresponding wave packet, 
but, in contrast to the wave packet, without dispersion. This group velocity is 
much smaller than the sound velocity. Finally, DB of type III are considered, 
being strongly localized and having an amplitude which is not small. The well- 
known Page and Sievers-Takeno modes are of this type, and numerically they have 
frequently been observed in a traveling state. It should be possible to provide an 
analytic estimate of their traveling speed by constructing an effective Hamiltonian, 
following the procedure of Section O and making use of Ansatz ()54|l . but this is yet 
to be done. The important feature which distinguishes DB of type III from those 
of type II is that their Peierls-Nabarro barrier is non-negligible. This barrier is 
estimated without explicitly computing the effective Hamiltonian, and the estimate 
is in reasonable agreement with the energy necessary to render a DB mobile in a 
numerical experiment. 
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